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Abstract 

In this paper, new schemes to synchronize linearly or nonlinearly coupled chaotic 
systems with an adaptive coupling strength are proposed. Unlike other adaptive schemes, 
which synchronize coupled chaotic systems to a special trajectory ( or an equilibrium 
point) of the uncoupled node by adding negative feedbacks adaptively; here, adaptive 
schemes for the coupling strength are used to synchronize coupled chaotic systems 
without knowing the synchronization trajectory. Moreover, in many applications, the 
state variables are not observable; instead, some functions of the states can be observed. 
How to synchronize coupled systems with the observed data is of great significance. In 
this paper, synchronization of nonlinearly coupled chaotic systems with an adaptive 
coupling strength is also discussed. The validity of those schemes are proved rigor- 
ously. Moreover, simulations show that by choosing some proper parameter a, the 
coupling strength obtained by adaptation could be much smaller. It means that chaos 
oscillators can be easily synchronized with a very weak coupling. 

Key words: Synchronization, Adaptive Coupling Strength, Nonlinearly Coupling, 
Time-varying Coupling Matrix. 



1 Introduction 



Recently, an increasing interest has been devoted to the study of complex networks (see 
Strogaze 2001 - Wang et.al 2002). Among them, synchronization of coupled complex net- 
works has received more attentions, because synchronization not only can explain many 
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natural phenomena (see Mirollo and Strogatz 1990), but also have many applications, such 
as image processing, secure communication (see Wei and Jia 2002 - Lu and Chen 2004) and 
so on. 



Generally, linearly coupled systems can be described as: 



&i(t) = f(xi(t),t) + cJ2 a ij T 



Xj(t) - Xi(t) 



1,2,- -.,JV (1) 



where x^t) = (xj(t), • • ■ , x™(t)) T e R n ; /(-, t) : R n x R + -> R n is continuous. The outer 

JV 

coupling matrix A satisfies: > 0, for i ^ j, and let an — — a ij\ while the inner 

coupling matrix T = diag^i, ■ ■ ■ , 7„) is positive definite, c denotes the coupling strength. 

If hm \\xi(t) — Xj(t)\\ = for all i,j = 1, 2, • • • , N, where || • || denotes some norm, 
then coupled systems dU) are said to be synchronized completely. In the following, we 
will investigate complete synchronization. Hitherto, many approaches and criteria to ensure 
synchronization have been derived (Lu et. al 2004 - Pikovsky, Rosenblum, and Kurths 
2001). 

Moreover, in practice, the state variables Xi(t) may be unobservable; instead, we can 
observe g(xi(t)), where g(-) is a monotone increasing function. How to synchronize cou- 
pled systems with the observed data is of great significance. In (Chen, Zhu 2007), following 
nonlinearly coupled systems 



Xi(i) = f(xi(t),t) + cJ2< 



g{xj{t)) -g{xi{t)) 



1,2,---,N (2) 



were proposed, where g(xj(t)) = (gi(xj(t)), ■ ■ ■ , g n (x^(t)) T and every is a nonlinear 
monotone increasing function. 

It is pointed out in Pikovsky, Rosenblum, and Kurths 2001 that if the coupling strength 
is larger than a critical value c*, then coupled systems <[TJ can be synchronized. In fact, c* 
depends not only on the coupling matrix but also depends on the dynamical behavior of 
function f(x(t),t) in the uncoupled system x(t) = f(x(t),t). Thus, the critical value c* 
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(suitable for all systems with various f(x(t),t)) is much larger than the coupling strength c 
needed for specified coupled systems (see Lu and Chen 2004). 

On the other hand, it is well known that a chaotic attractor typically has embedded 
within it an infinite number of unstable periodic orbits. In (Ott, Grebogi and Yorke 1990), 
it is shown that one can convert a chaotic attractor to any one of a large number of possible 
attracting time-periodic motions by making only small time-dependent perturbations of an 
available system parameter. On the other hand, if the attractor is not chaotic but is, say, 
periodic, then small parameter perturbations can only change the orbit slightly. 

Therefore, how to synchronize a large number of specified chaotic oscillators (for ex- 
ample, Lorenz oscillator or other oscillators) with a relatively small coupling strength is an 
interesting problem. 



For this purpose, we replace coupled systems ([T]) and © by 



ii(t) = f(xi(t),t) + c(t) J2 a ij T 



Xj(t) - Xi(t) 



g(xj(t)) - g(xi{t)) 



1,2,- ..,JV (3) 



l,2,---,iV (4) 



x i (t) = f(x i (t),t)+c(t)J2> 

with an adaptive coupling strength c(t). 

The adaptation of parameters is widely used in the signal processing and other research 
fields. Recently, the adaptive approach is used to synchronize master-salve systems or mu- 
tually coupled systems to a specified trajectory of the uncoupled node by adding negative 
feedbacks (see Chen and Zhou 2006 - Boccaletti 2006). For example, in (Chen and Zhou 
2006 and Zhou et.al 2006), the authors investigated how to synchronize coupled systems: 

±i(t) = f(x i (t),t) + h i (x 1 (t),x 2 (t),---,x N (t)) l<t<N (5) 



to a specified trajectory s(t) = f(s(t),t) by adding controls —di(t) 
adaptation rule di — ki\\xi(t) — 



Xi(t) - s(t) 



with 
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However, this adaptation rule does not suit complete synchronization. Because, the 
synchronization trajectory, which generally is not a trajectory of the uncoupled system, is 
unknown. 

In this paper, new schemes to synchronize chaotic oscillators with an adaptive coupling 
strength are proposed. It is revealed that a large number of chaotic oscillators (including 
Lorenz oscillators, Chen's oscillators, Rossler oscillators and Chua's circuits) can be syn- 
chronized with a very small coupling strength. It means that one can synchronize a large 
number of chaotic oscillators with very weak coupling. 

This paper is organized as follows: In section 2, some necessary definitions, lemmas, 
and hypotheses are given; in section 3, adaptive schemes for linearly and nonlinearly cou- 
pled system are given; in section 4 and section 5, some simulations are given to verify our 
theoretical results; while in section 6, an example is given to illustrate that: for coupled 
periodic systems the coupling strength should be very large; and the paper is concluded in 
section 7. 

2 Preliminaries 

In this section, some definitions, denotations and lemmas throughout the paper are pre- 
sented. 

Definition 1 The set S= {(xj , x%, ■ ■ ■ , xJf) T \xi — Xj] i, j — 1, 2, • • • , N} is called the 
synchronization manifold. 

Definition 2 Matrix A = (aij)^ =1 of order N is said to satisfy Condition Al, if 

N 

1. dij >0,i 7^ j, an = - £ a,ij,i = 1,2, •• -,7V 

2. Eigenvalues of A are all negative except an eigenvalue with multiplicity 1. 
Furthermore, if A e Al, and i j, then we say A e A2. 
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Definition 3 Suppose that A = diag{5\, 62, ■ ■ ■ , S n } is a diagonal matrix, and w > 0. 
/ : R n x R + — > R n is continuous, we say f e QUAD (A, w), if and only if 

{x - yf[f{x, t) - f(y, t)} - (x - yfA(x - y) < -w{x - yf(x - y) (6) 

holds for any x, y e R n . 

Remark 1 It is some control for the case (x — y) T [f(x, t) — f(y, t)] > (roughly speak- 
ing, f(x,t) — f(y,t) and x — y have same sign), \\f(x,t) — f(y,t)\ \ should be less than 
(x — y) T [A — zul n ](x — y). Instead, if (x — y) T [f(x, t) — f(y, t)) < 0, then inequality (0) 
is satisfied automatically. 

Lemma 1 If a matrix A^ x n £ Al. Then (see Wu 2005, Lu Chen 2006): 

1. [1, 1, • • • , 1] T is the right eigenvector of A. corresponding to eigenvalue with multi- 
plicity 1; 

2. The left eigenvector of A: £ = [£1, £ 2 , • ■ • , 6v] T £ R N corresponding to eigenvalue 
has following properties: it is non-zero and its multiplicity is 1; all & > 0, i = 1, 2, ■ ■ • , N. 
More precisely, A is irreducible if and only if all £3 > 0, i — 1, 2, • • • , N. 

N 

In the following, we always assume that A is irreducible and X) & = 1, & > 0, which 

i=l 

is called the normalized left eigenvector. It is clear that if A^ x n £ A2, then £j = -h for 

i = l,. ..,iV. 

3 Main Results 

In this section, we give the main results about how to design synchronization algorithms 
with an adaptive coupling strength for linearly and nonlinearly coupled systems. 
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3.1 Linearly coupled systems 



In this part, we discuss synchronization of linearly coupled systems with the coupling 
strength adaptively. 



3.1.1 Constant coupling matrix 

In this part, we investigate the case that the coupling matrix is time independent. 
Consider following linearly coupled systems with an adaptive coupling strength: 



the coupling matrix A = (a^) is not assumed to be symmetric. 

For any irreducible (asymmetric) matrix A^ x n G Al, Denote H = diag{£i, ■ ■ ■ , £at}, 
U — H — ££ T , where £ = ■ • ■ , £at) t G i?^ is the normalized left eigenvector corre- 
sponding to eigenvalue 0.. It is easy to check that —U G A2. 

Furthermore, denote X(t) = (xj(t), • ■ ■ , x T N {t)) T , F(X{t)) = (f(xi(t), t) T , • ■ • , 
f(x N (t),t) T ) T , A = I N ® A, A = A <g> T, and U = U ® I n , S = H (8 /„, where ® is the 
Kronecker product. Then, © can be rewritten in a compact form: 



The following Lemma gives a necessary and sufficient condition whether a coupled 
system is synchronized. It is clear that the left eigenvector £ plays a key role to prove 
synchronization. 

Lemma 2 X{t) G S if and only if 



Xi(t) = f(Xi(t),t) + c(t) Oij-r Xj(t) - 



(7) 



X(t) = F(X(t)) + c(t)AX(t) 



(8) 



X T (t)XJX(t) = 



(9) 
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In fact, for any X(t) = (xj (t), • • • , x%(t)) T , Y(t) = (yf(t), ■ ■ ■ , y T N {t)) T , we have (also 
see Wu 2005) 

1 N 

X T {t)\JY{t) = - E Uj(-r;ih - Xj(t)) T ( yi (t) - yj(t)) (10) 

Now, we propose a new scheme to synchronize linearly coupled systems with an adap- 
tive coupling strength and prove the following theorem. 

Theorem 1 Suppose A G Al is irreducible and f G QUAD(A,w). Then, following 
coupled systems with an adaptive coupling strength c(t): 

X(t)=F(X(t)) + c(t)AX(t), 
c (t) = -|X T (i)HAX(t) 

can finally achieve synchronization with a relatively small coupling strength, where c(0) = 
and a > 0. 

Proof: Pick a sufficiently large constant c > and define a Lyapunov function as: 

V(X(t)) = \x T (t)\JX(t) + -(c- c{t)f (12) 

Then, noting UA = EA, we have 

^ ( ^ (t)) = X(t) T U{F(X(t)) + c(t)AX(t)] + (c-c(t))X T (t)~AX(t) 

= X(t) T U[F(X(t)) - AX(t)\ + X T (t)VAX(t) + cX T (t)BAX(t) 

By the identity (flOb and the assumption / G Q£7 AD (A, zu), we have 

X r (t)U[F(X(t)) - AX(t)] 

N 

JV 

< £ ^•[^(0-^(0] r [^(0-^(*)] = -^ T (*)ux(*) 



Therefore, 

dV(t) 



< -wX{t) T \JX{t) + X T (t)UAX(t) + cX T (t)SAX(t) (13) 



dt 

Write Xj(t) = (x\(t), ■ ■ ■ , x 3 N {t)) T for j = 1, 2, • • • , n, we have 

n n 

X T (t)U(A + cA)X(t) = ^xJ(t)^C/%(t) + c^7 i xJ(t)HA£ i (t) (14) 

j=i i=i 

It can be seen that EA + A T E is a symmetric matrix with negative diagonal and row- 
sum zero. Now, let Vi, • ■ ■ , v N be the normalized eigenvectors of the matrix l(EA + A T E) 
with corresponding eigenvalues Ai = > A 2 > • • • > Ajv. Moreover, by x*(t) denote the 
projection of Xj(t) on the subspace L spanned by v 2 , • • ■ , Ujv- Then, 

n n 

c^xjit^jEAxjit) = -J2xJ(t) lj (EA + A T E)i 1 (t) 
i=i - i i 

n n 

= o E^* (tMZA + A T E)x*(t) < cX^ljZfmit) (15) 

Z 3=1 3=1 
n n 

xj (t) Sj Uxj (t) = £ xf (t) S 3 u %*j (*) ( 1 6 ) 

3=1 3=1 

Substituting into (fT3l . we have 

(]V(t) nn n T n T 

— A2 < -roI(fUI(f) + J] (tfaUxffl + cA 2 ^ 7j x* (t)x*(t) (17) 

Ci6 



It is clear that if c is sufficient large, then 



n n 

Y,xf{t)SjUxffl+c\ 2 Y,l3xf(t)£W) < (18) 

3=1 3=1 

In summary, we have ^ < -zuX(t) T XJX(t) < 0. 

Similarly, it can be seen that c(t) > and c(t) = if and only if X G S, so c(t) > for 
t > 0. 
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It is obvious that V — if and only if X G S. According to the well-known Lyapunov- 
LaSall type theorem for functional differential equations (see Kuang 1993), the trajectory 
of coupled systems, starting with arbitrary initial value, converges asymptotically to the 



X T UX = 0,c = c e R + }. Therefore, X T (t)UX(t) -> and c(t) -> c for some constant 
c - 

Theorem 1 is proved completely. 
Remark 2 Adaptive algorithm is often used in many research fields. Many authors use 



positive constants and s(t) is a special solution (or an equilibrium point) of the uncoupled 
system. Here, the synchronization state is unknown. Therefore, previous negative feedback 
adaptation methods are invalid. Theorem 1 provides a new adaptive algorithm, which suc- 
ceeds in the synchronization of complex networks with an adaptive coupling strength. 

3.1.2 Unknown constant coupling matrix 

In some physical coupled systems, the coupling matrix may be unknown, though we 
know A G Al. Can we design an adaptive algorithm to synchronize coupled systems with 
an adaptive coupling strength for a unknown coupling matrix? In this subsection, we will 
give an affirmative answer by proving the following theorem. 

Theorem 2 Suppose that the unknown coupling matrix A G Al is irreducible and f G 
QU AD (A, w). Then, following coupled systems with an adaptive coupling strength c(t): 



largest invariant set Hi contained in H 2 = {V(t) = 0} as t — > +oo, where Hi = {[X, c] 



T . 



di Xi(t) — s(t) as the negative feedback adaptation, where di = ki\\xi(t) — s(t) \\ 




< 



X(t)=F(X(t)) + c(t)AX(t), 
c(t) = -f X T (t)kX(t) 



(19) 



can achieve synchronization with a relatively small coupling strength, where A — A <g> 7, 
A G A2 is any irreducible matrix, c(0) = and a > 0. 
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Proof: Define a slightly different Lyapunov function 

V^X®) = l -X T (t)VX(t) + l(c- c(t)) 2 (20) 
1 a 

where rj > is a constant. Then, noting UA — EA, we have 

dV l {X{t)) vttTTTfp/vm , , , _ u „ vT 



dt 



= X(ty V[F(X(t)) + c(t)AX(t)\ + 77(0 - c(t))X 7 (t)AX(t) 
= X(t) T U[F(X(t)) - AX(t)\ + X T (t)VAX(t) 
+c(t)X T (t)SAX(t) + r](c - c{t))X T {t)AX{t) 
< -mX(t) T XJX(t) + X T (t) (UA + cr]A^X(t) 



+c{t)X T {t) (sA-r]A^X(t) 



Now, pick a sufficiently large constant c and a sufficiently small 77 such that 

UA + crjA < 
{HA} S - 77A < 



(21) 



Therefore, dVl ^f^ < —zuX(t) T XJX(t) < 0. By similar arguments used in Theorem 1, we 
can prove Theorem 2. 



3.1.3 Time-varying coupling matrix 

In practice, the coupling matrix is time-dependent, which means that the coupling matrix 
changes along with time. Therefore, it is natural to investigate following linearly coupled 
systems with a time-varying coupling matrix: 

X(t) = F(X(t)) + c(t)A(t)X(t) (22) 



The following theorem for time-varying coupling can be proved with the same argu- 
ments used in the proof of Theorem 1 . 
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Theorem 3 Suppose that f G QUAD(A,w), the coupling matrix A(t) G Al is irre- 
ducible and has the same left eigenvector £ corresponding to eigenvalue for all t. If the 
largest non-zero eigenvalue of the matrix EA(t) + A T (t)E satisfies \2(t) < A < 0/or all t. 
Then, following coupled systems with an adaptive coupling strength c(t): 

X(t)=F(X(t)) + c(t)A(t)X(t), 
c(t) = -fX r (t)EA(t)X(t) 

can finally achieve synchronization with a relatively small coupling strength, where c(0) = 
and a > 0. 



Remark 3 The condition that A(t) G Al w irreducible and has the same left eigenvector 
£ corresponding to eigenvalue Ofor all t looks quite strong. However, it just suits the node- 
balanced coupling networks discussed in (Belykh, V. et.al 2006), where the coupling matrix 
A(t) is assumed to be row-sum zero as well as column-sum zero for each t. In this case, 
A{t) G Al has the same left eigenvector [1, • • • , 1] T corresponding to eigenvalue for all 
t. Therefore, Theorem 3 applies to the node-balanced coupling networks. 

As a direct consequence of Theorem 3, we can obtain the following simple adaptive 
scheme, if every A(t) is symmetric. 

Corollary 1 Suppose that f G QUAD(A, w), the coupling matrix A(t) G A2 is irre- 
ducible. If the largest non-zero eigenvalue of the matrix A(t) satisfies \2(t) < A < 0, for 
all t. Then, following coupled systems with an adaptive coupling strength c(t): 

' X(t)=F(X(t)) + c(t)A(t)X(t) : 
c(t) = -fX T (t)A(t)X(t) 

can finally achieve synchronization with a relatively small coupling strength, where c(0) = 
and a > 0. 
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Remark 4 In Theorem 3 and Corollary 1, condition \2{t) < A < plays a key role. 
However, calculating X 2 (t)for all t is impossible numerically. If all A(t) can be dominated 
by a constant matrix A, then we can synchronize chaotic oscillators directly using rtiij) . only 
replacing A in the adaptive algorithm by A. Following corollary explains it. 

Corollary 2 Suppose that f e QUAD(A, w), the coupling matrix A{t) 6 A2 is irre- 
ducible. If there exists a constant matrix A = (a^ ) e A2 such that > Oy(t), i ^ j hold 
for all t. Then, following coupled systems with an adaptive coupling strength c(t) 

X(t)=F(X(t)) + c(t)A(t)X(t), 
c(t) = -f X T {t)AX{t) 

can finally achieve synchronization with a relatively small coupling strength, where c(0) = 
and a > 0. 

Proof In this case, & = ^, which implies EA(t) = jjA(t), and 

N 

X T {t)A(t)X(t) = £ %•(*) - Xj {t)) T {xi(t) - Xj (t)) 

N 

< ]T a ij {x i {t)-x j {t)f{x i (t)-x j {t)) = X T {t)AX{t) (26) 

Using the Lyapunov function 

V(X(t)) = l -X T {t)\JX(t) + -(c- c{t)f (27) 
2 a 

we have 

< X{t) T U[F{X{t))- AX{t)]+X T {t)U{A + cNA)X{t) 
By the same arguments used in the proof of Theorem 1 . we can obtain Corollary 2. 
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3.2 Nonlinearly coupled systems 

We consider following nonlinearly coupled systems with an adaptive coupling strength: 

(28) 

-i 

where g( Xi (t)) = [gi(x}(t)), ■ ■ ■ , g n (x? (t))] T . 



ii(t) = f(xi(t),t) + c(t) a ij g( x j(t)) - g( x i(t)) 



Denote G(X(t)) = (g T (x!(t)), ■ ■ ■ ,g T (x N (t))) T and A = A <g> I n , then we have 

^ = F(X(t)) + c(t)AG(X(t)) (29) 

Theorem 4 Suppose A e A2 m irreducible, f G QC/AD(A, ca), gfc( ^If W > /3/or 
Vu ^ v,k = 1, ■ ■ ■ ,n, where (3 is a positive constant. Then, following nonlinearly coupled 
systems with an adaptive coupling strength c(t): 

X{t) = F{X(t)) + c{t)AG(X(t)), 
c(t) = -%X T (t)AG(X(t)) 

will be synchronized with a relatively small coupling strength, where c(0) = and a > 0. 

Proof: Because A e A2, then ^ = and £L4 = A Using the same Lyapunov 
function as in the Theorem 1 , we have 

= X(t) T V[F(X(t)) + c(t)AG(X(t))] + ^-^))X T (t)VAG(X(t)) 

= X(t) T U[F(X(t)) - AX(t)} + X T (t)VAX(t) + cNX T (t)XJAG(X(t)) 

Furthermore, 

n n 1 

X T (t)UAG(X(t)) = x;sf = E ^(0^(^(0) 

1=1 i=l iv 

= - E 4 E - ^(*)] r [ftW(*)) - mtfm 

i=l JV j>fc 

< -pE 4 E - K (*) - *?(*)] 

1=1 iv j>fc 

n i n 

= /?E A? 5 *(*) T ^i(*) = /5E^(*) T ^<(*) 

1=1 iv 1=1 
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where <7i(£i(f)) = (gi(x}(t)), ■ ■ ■ , g^xf (t))) T , fori = 1,2, • • • ,n. Therefore, 
dV{t) 



] < -wX{t) T XJX{t) + Y,xj(t)[6jU + c(3NUA\ Xj {t) (31) 
at j=1 

By similar arguments used in Theorem 1, we complete the proof of Theorem 4. 

Similar to linearly coupled systems, for nonlinearly coupled systems with a time- varying 
coupling matrix: 

dX(t) 
dt 

we have 



F(X(t)) + c(t)A(t)G(X(t)) (32) 



Corollary 3 Suppose A{t) e A2 is irreducible, f G QUAD (A, w), ih^hMhM > pf or 
Vm 7^ v,k = 1, • • • , n, where f3 is a positive constant. If the largest non-zero eigenvalue of 
A(t) satisfies A 2 (t) < A < 0, for all t. Then, following nonlinearly coupled systems with an 
adaptive coupling strength 

X(t) = F(X(t)) + c(t)A(t)G(X(t)), 
tit) = -%X T (t)A(t)G(X(t)) 

will be synchronized with a relatively small coupling strength, where c(0) = and a > 0. 



4 Numerical Simulation 

To validate the effectiveness of the proposed synchronization algorithms with an adaptive 
coupling strength, in the following, we couple 100 chaotic oscillators by following linearly 
coupled systems with adaptive strength 

X(t)=F(X(t)) + c(t)AX(t), 
c(t) = -^X T (t)EAX(t) 
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and nonlinearly coupled systems with adaptive strength 



< 



X(t) = F(X(t)) + c(t)AG(X(t)), 
c(t) = -ZX T (t)BG(X(t)) 



(35) 



where X(t) = (xj(t), ■ ■ ■ , xf 00 (t)f, F(X(t)) = (f( Xl (t), t) T , ■ ■ ■ , f(x wo (t), t) T ) T , A = 
A®I 3 , A e i? 10 ^ 100 . b = B®h, and B = \{A+A T ), G(X(t)) = (<?(xi(£)), g(x 2 (t)), 
^(^ioo(^))) T ,^(^W) = (x}(t) +tanh(x\(t)),xj(t)+tanh(xj(t)),xl(t)+tanh(xl(t))) T , 



i = 1, • • • , 100. The prototypes of /(•) are Chua's circuit, Chen's oscillator, Lorenz's oscil- 
lator and Rossler's oscillator, respectively. For the coupling matrix A, first, we construct a 
coupling matrix of a small-world network generated by the method proposed in (Watts et.al 
1998). Then, replace every non-zero element a^-, % ^ j with a positive random scalar, which 
is equally distributed in [0, 1], and choose the diagonal elements to ensure A e Al. 

We also couple 100 chaotic oscillators by following linearly coupled systems with an 
adaptive strength 



In all simulations, we use the quantity E(t) = J J2 11^(0 — ^i(^)|| 2 /99 as a measure 



X(t)=F(X(t)) + c(t)AX(t), 
c(t) = -ZX T (t)AX(t) 



(36) 



where A = A <g> I 3 , A e R 
matrix respectively. 



100x100 



are globally connected matrix or randomly connected 




of synchronization error. 



4.1 Chua's Circuits 



The uncoupled equation is : 



dx 
dt 
dy 
~dt 
dz 
dt 



x — y + z 



m[y - h(x)] 



-ny 
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where h(x) = Ix — j^[\x + 1| — \x — 1|], fh = 9 and ri = 14=. Figure 1. (a) and (b) show 
the dynamics of c(t) and E(t) for linearly coupled systems d34l) : while (c) and (d) show the 
dynamics of c(t) and E(t) for nonlinearly coupled systems d35t . 




c=0.031935 



10 20 
(a) Time t 



30 




10 20 
(b) Time t 




c=0.032586 



10 20 
(c) Time t 



30 




10 20 
(d) Time t 



Figure 1. The dynamics of c(t) and E{t) for 100 linearly and nonlinearly 
coupled chua's circuits with an adaptive coupling strength 

Figure 2. (a) and (b) show the dynamics of c(t) and E(t) for linearly coupling systems 
d36t when A is the globally coupled matrix; while (c) and (d) show the dynamics of c(t) 
and E(t) for linearly coupled systems d36t when A E A2 is a random matrix. 
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10 20 30 10 20 

(c) Time t (d) Time t 



Figure 2. The dynamics of c(t) and E(t) for d36l) with globally and randomly connected A 

4.2 Chen's Oscillator 

The uncoupled equation is : 

x = a(y — x) 
< y = (c — a)x — xz + q/ 
i = xy — bz 

< 

where a = 35, b = 3 and c = 28. Figure 3. (a) and (b) show the dynamics of c(t) and E{t) 
for linearly coupled systems d34t : while (c) and (d) show the dynamics of c(t) and E(t) for 
nonlinearly coupled systems (l35l) . 
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Figure 3. The dynamics of c(t) and E{t) for 100 linearly and nonlinearly 
coupled Chen's chaotic oscillators with an adaptive coupling strength 

Figure 4. (a) and (b) show the dynamics of c(i) and E(t) for linearly coupling systems 
d36t when A is the globally coupled matrix; while (c) and (d) show the dynamics of c(t) 
and E{t) for linearly coupled systems d3*6l ) when A G A2 is a random matrix. 
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Figure 4. The dynamics of c(t) and E(t) for d36l) with globally and randomly connected A 

4.3 Lorenz's Oscillator 

The uncoupled equation is : 

xi = (3(x 2 - xx) 
< X2 = axi — X1X3 — £2 
x 3 = X1X2 — bxs 

where /3 = 10, a = 28, and b = |. Figure 5. (a) and (b) show the dynamics of c(t) and E(t) 
for linearly coupled systems (l34b : while (c) and (d) show the dynamics of c(t) and for 
nonlinearly coupled systems (l35l) . 
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Figure 5. The dynamics of c(t) and E{t) for 100 linearly and nonlinearly 
coupled Lorenz's chaotic oscillators with an adaptive coupling strength 



Figure 6. (a) and (b) show the dynamics of c{t) and E{t) for linearly coupling systems 
d36t when A is the globally coupled matrix; while (c) and (d) show the dynamics of c(t) 
and E{t) for linearly coupled systems ( BBT ) when A 6 A2 is a random matrix. 
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Figure 6. The dynamics of c(t) and E(t) for d36l) with globally and randomly connected A 



4.4 Rossler's Oscillator 



The uncoupled equation is : 

X\ = — %2 — 
' X 2 = Xi + 0.2X2 

x 3 = 0.2 + x 3 (xi — ji) 

where /i = 5.7. Figure 7. (a) and (b) show the dynamics of c(t) and E{t) for linearly 
coupled systems (|34b : while (c) and (d) show the dynamics of c(t) and E(t) for nonlinearly 
coupled systems (l35l) . 
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Figure 7. The dynamics of c(t) and E(t) for 100 linearly and nonlinearly 
coupled Rossler's chaotic oscillators with an adaptive coupling strength 

Figure 8. (a) and (b) show the dynamics of c(i) and E(t) for linearly coupling systems 
d36t when A is the globally coupled matrix; while (c) and (d) show the dynamics of c(t) 
and E{t) for linearly coupled systems ( BBT ) when A G A2 is a random matrix. 
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Figure 8. The dynamics of c(t) and E{t) for d36l) with globally and randomly connected A 



5 Numerical Example 2 



In this section, we give some simulations for coupled systems with a time-varying cou- 
pling matrix. As a special case, it includes node-balanced coupled systems. 



We pick 



A(t) 



Pi 



P3 J 



-5 — sin t — cos t 3 + sin t 2 + cos t 

2 + cost — 5 — sin t — cost 3 + sint 

3 + sint 2 + cost — 5 — sin t — cost 



(37) 
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as the time-varying coupling matrix, where pi,p 2 ,P3 are positive constants. 

It is clear that for all t, A(t) has the same left eigenvector £ = (^-, ^-) T corre- 
sponding to eigenvalue 0. In particular, if pi = p 2 = ps, then the coupling matrix A(t) is 
node-balanced for each t. 

Moreover, it is easy to check that 

\ 2 (t) = -(5 + sint + cost) < -(5 - V2) < 



is the largest non-zero eigenvalue of the matrix 



,,,Tr, 5 + sint + cost 
EA(t) + A(t) T E = 



-2 1 1 ^ 
1 -2 1 
1 1 -2 j 



We couple three Chua's circuits by 



X(t) = F(X(t)) + c(t)A(t)X(t), 
c(t) = -|X T (t)HA(t)X(t) 



(38) 



We also use E(t) 
synchronization. 



\x 2 (t) — ^i(t)|| 2 + ||^3(0 — x i{t)\\ 2 )/2 to measure the error of 



Figure 9. shows the dynamics of the coupling strength c(t) and error for different 
Pi,P2,P3, Where sub-figures (a) and (b) denote the dynamics of c(t) and E(t) for p 1 = 
p 2 — 1, Pz — 2; and sub-figures (c) and (d) denote the dynamics of c(t) and E(t) for 

Pl=P2=P3 = 1- 
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Figure 9. The dynamics of c(i) and E(t) for the coupled system (l38l) . 

6 Conclusions 

In this paper, we propose new algorithms to synchronize linearly or nonlinearly coupled 
systems with an adaptive coupling strength. Unlike those adaptive algorithms existing in 
the literature, where coupled systems are synchronized to a special trajectory s{t) or an 
equilibrium of the uncoupled system by adding a negative feedback controller; in this paper, 
we synchronize linearly and nonlinearly coupled systems with an adaptive coupling strength 
without knowing the synchronization trajectory. By adapting the coupling strength, we 
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reveal that a large scale of chaotic oscillators can be synchronized even with a very small 
coupling strength. It indicates that chaotic oscillators are very easy to be synchronized. 
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